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suriace. l nis leads to an integro-ditierentiai operator m tne equation tor tne 
radial distribution of the mean magnetic field strength, Q{r) in the disc plane 
at a distance r from its centre; an asymptotic form of its solution at large 
distances from the dynamo active region is obtained. Numerical solutions of 
the integro-differential equation confirm that the non-local effects act similarly 
to an enhanced magnetic diffusion. This leads to a wider radial distribution 
of the eigensolution and faster propagation of magnetic fronts, compared to 
solutions with the radial surface field neglected. Another result of non-local 
effects is a slowly decaying algebraic tail of the eigenfunctions outside the 
dynamo active region, Q{r) ~ r~^, which is shown to persist in nonlinear 
solutions where a-quenching is included. The non-local nature of the solutions 
can affect the radial profile of the regular magnetic field in spiral galaxies and 
accretion discs at large distances from the centre. 

KEY WORDS: Mean-field dynamos, thin-disc asymptotics, boundary condi- 
tions, galactic magnetic fields, accretion discs. 

1 Introduction 

Thin-disc asymptotics, applicable to spiral galaxies and accretion discs, have been a 
useful tool in the studies of the origin of large-scale magnetic fields in these objects. 
The small parameter naturally arising in a thin disc is its aspect ratio 

A = ho/Ro < 1 , (1) 

where ho and Rq arc the characteristic half-thickness and radius of the disc respec- 
tively. The height of the disc surface and the radius of its edge, understood as 
positions where suitable boundary conditions for the disc's magnetic field can be 
reasonably applied, are h — O{ho) and R — 0{Ro)- In a domain that has no sharp 
boundaries (e.g., a plasma layer in hydrostatic equilibrium along the vertical direc- 
tion and centrifugal equilibrium along radius), ho and Ro can be identified with the 
scale height and the radial scale length, respectively. 
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we aaopt cyimaricai poiar cooramates [r,(p,z), m wnicn z measures aistance 
parallel to the rotation axis and assume that our system is axisymmetric, indepen- 
dent of the azimuthal angle 0. We assume that the mean motion is dominated by 
differential rotation and is approximated by the velocity V = (0, V, 0). We con- 
sider an axisymmetric magnetic field whose poloidal components can be expressed 
in terms of the azimuthal component A of the vector potential: 

B.r-|i.5.i4M)V (2) 



dz^ ^ r dr 

The magnetic field satisfies the mean-field dynamo equation 

— = V X (V X B + aB - /3V X B) (3) 

(see Moffatt, 1978; Krause and Radler, 1980), in which a and P are turbulent 
transport coefficients responsible for the a-effect and turbulent magnetic diffusion, 
respectively. We will consider the linear problem with V, a and (3 assumed given and 
independent of time t, for which solutions may be sought proportional to exp(rt), 
where the constant F is the growth rate. 

The dynamo equation (jS)) has been solved with a non-local boundary condition, 
possible for an infinite slab, using the Bullard-Gellman (1954) formalism (Raedler 
and Wiedemann, 1990). Traditionally, however, the boundary conditions used in 
conjuction with (jS)) are the so-called vacuum boundary conditions. If there are no 
electric currents outside the disc, then V x B = there, and the magnetic field may 
be expressed as the gradient of a potential in the form B = — V$. Thus the axial 
symmetry implies that the azimuthal magnetic field vanishes, 

B = outside the disc (4) 

(see, e.g., Zeldovich et ai, 1983, p. 151), while the remaining meridional magnetic 
field is potential = since V-B = 0. Our solution of the potential problem relies 
on two approximations. Firstly, we assume that the disc is so thin that its surface 
is essentially at z = 0. In this way we may obtain, using the Green's function. 



a solution oi tne i\eumann proDiem lor ^?[r,z) m tiie upper nan plane z > u m 
which the value of ^^(r, 0+) = — 9$/9z(r, 0+) is specified. This Green's function 
solution may by used to calculate Br{r, 0+) = —d^/dr{r, 0+). In terms of the vector 
potential A for the radial and axial magnetic field components Br = —dA/dz and 
Bz = {l/r)d{rA)/dr, the resulting expression for Br in terms of Bz, both at 2; = 0+, 
may be written in the form 

dA 



dz 



(r,0+)-£{A(r,0+)} = 0^ 



where the operator 



/:{A(r,0+)} 



involves the Green's function 



d_ 

dr 



G{ry)-^[r'A{r\0^)) 



dr' 



de 



i/r^ — 2rr' cos 9 + r'^ 



(5) 



(6) 



(7) 



Secondly, we assume that the magnetic field is localised above the disc and is neg- 
ligible beyond r = R. That means that the upper limit 00 in the integral © is 
interchangeable with R. Then (jSj) is to be regarded as an equation that determines 
Br as a nonlocal integral function of Bz on the disc surface z = 0+ for vacuum 
boundary conditions. Indeed for the quadrupolar magnetic fields, which concern us 
in this paper, we also have Bz{r, 0) = for r > R. For that case, no approximations 
are involved in interchanging 00 with R. 

In a previous paper, Priklonsky et al. (2000, hereafter referred to as Paper I) 
adopted the slightly different formulation 



C{Air,0^)}=l j^W{ry) A ^ (rM(r', 0^)) 



dr' 



of the integral operator based on the solution of (V^ — r^^)A = in the vacuum 
region outside the disc, where W{r,r') will be given below, see It is readily 

established by integration by parts that the operators © and (jH)) are equivalent 
provided that the two kernels G{r,r') and W{r,r') are related by 



d 



1 d 



r^[G(r,r')] = - - — [l^(r, r')] . 



(9) 



it IS trie latter lormuiation (j^ tnat we employ to oDtam our numerical results. 
The numerical implementation involves a regularisation procedure to allow for the 
singularity of the kernel, W ~ — Tr^^v^rr'ln \r — r'\ at r = r' . 

Inside thin discs, \z\ < h{r) on < r < R, where h may be a function of r, WKB 
type solutions of (jHI) can be constructed in the form 



(Soward, 1978, 1992a,b; Ruzmaikin et ai, 1985, 1988). The local z-structure on the 
short length scale h of the disc thickness, is determined by (6, a) normalised at our 
convenience; here z is measured in units of h so that the boundary is located at 
z = ±1. The corresponding time unit used to measure the inverse growth rate l/F 
is /io//5- The modulation on the long radial length scale i^o, adopted as our unit 
of r, is quantified by the amplitude Q{r). In our thin disc (A <ti 1) the solution 
Q{r/\^) is, in fact, usually modulated on a radial scale X^Rq intermediate between 
i?o and /lo, where s is a constant satisfying < s < 1. The scale X'^Rq is large 
compared to the disc width, h = 0{XRq), but small compared to the disc radius, 
R, a quantity of order Rq. The actual value taken by s depends on the nature of 
the boundary conditions. It is s = 1/3 for our non-local boundary condition © but 
may take other values if different boundary conditions are adopted; e.g., the value 
s = 1/2 found by Ruzmaikin et al. (1988) for a 'local' boundary condition as we 
explain below in relation to ()14|1 . 

We note that generally the growth rate F of the dominant modes in a thin disc 
is real, when ±a > for ±z > 0. Raedler and Wiedemann (1990) considered non- 
local boundary conditions for an infinite slab, also finding that non-oscillatory modes 
are the most easily excited. The asymptotic expansion of the form ()1()|1 applies to 
modes of quadrupolar parity that dominate in galactic discs; dipolar modes, perhaps 
prevailing in accretion discs (e.g.. Moss and Shukurov, 2003), have been discussed 
by Soward (1978, 1992a,b). 




with 




(10) 



5 



iiie lowest-oraer approximation m a yieias a Dounaary value proDiem lor [o,a) 
that only involves derivatives in z and so is local in r. The local problem yields the 
local eigenfunction (6, a) together with its eigenvalue, namely the 'local growth rate' 
7(r) (see Paper I for details). Relative to local coordinates, in which the units of 
length for z and r are h and Ro respectively {z z/h and r r/Ro), the local 
eigenfunction determines a linear relation between da/dz{l;r) and a(l;r) at the 
disc surface. We use dimensionless coordinates z ^ z/h and r — > t/Rq, keeping 
the notation z and r for them, unless stated otherwise. Used in conjunction with 
the boundary condition (jSJ, the linear relation leads at next order to an amplitude 
equation governing Q{r). Written in the form derived in Paper I, it is 

[r - 7(r)]g(r) = Xr]{r)C{q{r)}, (11) 

where the factor A reflects the fact that the radial unit of distance Rq is a factor 1/A 
larger than the axial unit h, the integral operator C is defined by (jHl) with kernel 
W{r,r') or © with kernel G{r,r') and 



q{r) = g(r)a(l; r) , r/(r) = p^^^ ' ^ 



(12) 



b{z; r) 
a{z; r) 

Here X is the eigenvector of the lowest-order boundary value problem, the asterisk 
denotes the eigenvector of its adjoint problem, while 

(X,X,) = / X-X,dz. 
Jo 

The solution of (fTT|) subject to the boundary conditions 

g(0) = and g ^ as r | cxo (13) 

provides yet another eigenvalue problem, for which the eigenvalue is the growth rate 
F and the eigenfunction is q{r) which determines the radial r-modulation Q. 

Soward (1978, 1992a,b) developed a systematic asymptotic theory for the marginally 
stable dynamo eigensolutions in a thin slab that includes the non-local coupling be- 
tween its distant parts via the surrounding vacuum, as described by the integral 



operator m ifiiji out m simpiinea locai uartesian ratner tnan axisymmetric lorm. 
Though this theory was originally formulated in terms of Fourier transforms, real 
space formulation is given in as obtained in Soward (2003). He pointed out, 
in particular, that the non-local coupling produces slowly decaying (algebraic) tails 
in the magnetic field radial distribution far away from the localization region of an 
eigenf unction. The corresponding theory with a full account for the cylindrical ge- 
ometry and with allowance for growing (rather than steady) solutions was proposed 
in Paper I. 

The integral term in (fTT|) was neglected by Ruzmaikin et al. (1988), who were 
interested mainly in solutions within the dynamo active region, where the non-local 
external potential coupling may not be very important, so the integral term in ^ 
was neglected. Instead these authors included local internal coupling by diffusion, 
which leads to the radial dynamo equation 



similar to ()11|). but with the integral term replaced by the diffusion operator. The 
value of the exponent s in ()1()|1 is sensitive to this difference. If (fT^ is used, then 
s = 1/2 (Ruzmaikin et al, 1985, 1988), whereas non-local asymptotics have s = 1/3 
(Soward, 1978). 

In this paper we discuss numerical solutions of the non-local radial equation (fTT|l . 
both in the kinematic regime and in a nonlinear regime with heuristically chosen 
nonlinearity. We compare solutions of the non-local equation (fTT|) with those of the 
local equation (fT^ . where the integral term is replaced by the diffusion operator, 
and confirm that the non-local coupling produces effects similar to enhanced radial 
diffusion. We also demonstrate that the algebraic tail of magnetic field outside the 
dynamo active region persists in nonlinear solutions. 

Before continuing it must be emphasised that though it was convenient in ()12|) 
to measure the z-distance in units of h for the internal disc magnetic field, this is 
quite inappropiate in the vacuum region outside. So for the remainder of our paper 




(14) 
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Dotn z ana r wui oe measurea m units oi rto? tnougn we continue to write a\i^r) 
for the value of a on the disc surface. 



2 The Integral Kernel 

In this section we develop various forms for the equivalent integral operators (jHl) and 
(jHI) employed in which may be utilised effectively in our numerical methods as 
well as aid our understanding of the asymptotic results. 



2.1 Properties of G(r, r') and VF(r, r') 



We begin by noting that the kernels G{r,r') and W{r,r') in (jH)) and (jH)) have the 
integral representations 



G{r, r') 



Jo{kr) Jo{kr') dk , 



W{r, r') 



rr I Jiikr) Jiikr') dk (15) 
'o 



in terms of the Bessel functions Jq{x) and Ji{x) of the first kind. The equivalence 
of the former for G{r,r') with ((Zj) follows from Gradshteyn and Ryzhik's (2000), 
(3.674.1) and (6.512.1), and Abramowitz and Stegun's (1965), (17.3.9); henceforth 
referred to as G&R and A&S respectively. Further it is readily established using 
the elementary properties of the Bessel functions and their derivatives that the two 
kernels G{r,r') and W{r,r') are related by 

Significantly the kernels have the alternative representations 



G{r,r') = —K{m) 



Tcr 



2r 



W{r,r') = — m 



TT 



dK 

K{m) + 2(m - 1)-— (m) 
dm 



for r' < r 



(16) 



in terms of the complete elliptic integral of the first kind 

f.7r/2 _7/i / f\ 2 



K{m) 



de 



and 



(17) 



'0 V 1 — msin^ 9 

(see G&R (6.512.1), A&S (15.2.7), (17.3.9)). Direct verification that G'(r,r') and 
W{r, r') defined by ()16|) satisfy the relations follows from the identity A{d/ dm) [m(l- 
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m)aj\ / am\ = a [see Lr&cn. (^8.iz4.ijj. ine lorms lor r > r are oDtamea using tne 
symmetry relations 



G(r,r') = G(r',r) 



iy(r,r') = W{r',r) . 



The asymptotic representations of the kernels are readily determined by docu- 
mented properties of the complete elliptic integral. In particular use of the formula 
2K{m) ^ ln[16/(l - m)] valid when 1 - m< 1 (see A&S (17.3.26)) shows that the 
kernels G{r,r') and W{r,r') have the expansions 



G(r, r') 
W{r, r') 



1 


In 


f 8r 






7rr 






r 


In 


f 8r 






TT 




, ^ ry 



for < r — r' ^ r . 



(19) 



Likewise use of the formula K{m) ^ (7r/2)(l + m/4) valid when m ^ 1 (see A&S 
(17.3.11)) shows that the kernel G{r,r') has the expansion 

.'2 



G{r, r') 



1 r 

^ — ^ 



for r' <^ r . 



(20) 



Together these results determine the important asymptotic forms 

1 



drdr 



■j{r, r') 



-Kr(r — r'Y 

3r' 

2^ 



for < r — r' <^ r , 
for r' <^ r . 



(21) 



2.2 Asymptotic properties of localised solutions 

In this section we discuss asymptotic forms of the solutions far from the dynamo 
active region. In particular, we clarify the nature of the slow, algebraic (rather than 
exponential) decay of the solution with cylindrical radius at large radii. 

When g(r) is localised about some point on a length scale e, where e -C r^, 
we integrate (jH)) by parts to obtain [noting that A{r, 1) = q{r) exp (Tt)] 



£ r'q{r')^\G{ry)]dr' 



(22) 
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2^™^-^ — {r,rm) for |r-r^|>e, (23) 
oror-r 



m 



where the slashed integral sign indicates principal part [needed here because dG/dr' ^ 
(r — r')~^] and 

/•oo 

Q = q{r) dr . (24) 
Jo 

This result may be used in conjunction with (fTT|) to obtain the far field behaviours 



q{r)^-- ^rmTT^ — (^,^m) for r-'^m|>e (25) 

i — 7(rj orOTn 



m 



of the modulation amplitude and the radial magnetic field 

1 dA d^G 
Br{r, 1) = -J-Q^if^^ 1) ~ exp(n)Qr,n ^^^^ (r,r^) for |r - r^] > e (26) 

on the disc surface, as defined by (jH))- Furthermore from ()2ip we may employ the 
limiting forms 

1 



for rm. > k — Tml > e , 



c^^G , , I 7r(r - r^)2 

The implication of the r-dependencies in and (j^7|) is that the response of the 
magnetic field in vacuum outside the disc corresponds to a line dipole source when 
rm \r — ^ e but corresponds to a quadrupole source when r ^ rm- 
To illustrate the idea we consider the 'local growth rate' 

l(,r) = 7^ - 1 |7^| (r^ - r)^ on < r < oo , (28) 

so that the infiuence of 7(r) is dominated by its behaviour near r^, where 7^ and 
7^ are constants. Then the amplitude q{r) governed by (fTT| satisfies 

[(r - 7n.) + i WL\ (rm - tY] q{r) = Ar^^£{g(r)} , (29) 
where rjm = vi^m)- Upon setting 

r = r^ + ezu, T = 7^ + ie^ 6 with e = {2Xr]J\i:^\y/\ (30) 
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tne integral equation ifz^jji may De approximatea lor w = U[i) witn tne neip oi 
and (ES) to yield 



{6 + zu')q{zu) ^ - r dw' , (31) 

TT J _^ w — w aw 

where q{w) = q{r). Here, since e <^ 1, the integration over < r < oo is equivalent 
asymptotically to — oo < < oo on the basis that g ^ as \w\ = \r — rm\l^ — ^ oo. 

The lowest eigenvalue 5 of the eigenvalue problem ()31|) given by Soward (2003) 
is 5 ~ —1.0188, namely the first zero of the derivative of Airy's function, Ai'((5) = 0. 
For that value we have the eigensolution 

^(^) — — t-\ / Ai(fc + (5) cos fctu (i/c , (32) 
7rAi(()) Jo 

where Ai(5) ^ 0.53565 and 

/oo 
q{w) dw . (33) 
-oo 

The large vj behaviour of q is 

(34) 

This result is consistent with fOKjl and (fTfjl provided that ^ |r — r^l ^ e. 



3 Numerical results 

In this section we discuss numerical solutions of the non-local radial equation (jllll 
and assess the importance and implications of the non-local effects by comparing 
the solutions with those of its local 'approximation', (fT^ . Our numerical solution 
of (jllll employs the integral kernel iy(r, r') defined by ()lfi|l . To that end K(m) 
can be computed iteratively using the procedure of the arithmetic-geometric mean 
described in Sect. 17.6 of A&S, which displays phenomenal convergence. We also 
need to compute dK/dm\ for this, we use a modification of the above procedure 
described in Appendix A. Numerical solutions of the linear equation (fTTjl and its 
nonlinear version, which we introduce in Section ESI involve a further regularization 
procedure described in Appendix B to contend with the singularity of the kernel. 
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iiie local aynamo proDiem at a nxea r, tiiat yieias tiie locai growtii rate 'j[r) ana 
the local eigenfunction r), a(z; r)), will not be solved here. Instead a suitable 
local growth rate, 7(r), is prescribed and we take 



a(l;r) = 1 



r]{r) = r]{r) = 1 



so that q{r) = Q{r) . 



(35) 



We adopt the boundary conditions 



g(0) = Q{R) = 0, 



(36) 



where R is the disc radius. Since our nonlocal integral equation (jllll is derived for 
a disc of infinite radius, this may lead to an artificial boundary layer a.t r = R if R 
is not large enough. Thus it is important that our magnetic field is localised and 
negligible for r > R. In the following subsections we consider examples of the cases 
in which 7(r) is maximised on the symmetry axis = and also at a finite distance 
> from it. 

3.1 Linear solutions for the case = 

The two leading eigenfunctions for the nonlocal (fTTj) and local (fT^ eigenvalue prob- 
lems are shown in Fig. Q for the parabolic profile 



with 7m = 10, = 0, ro = 1 and R = 2. The non-local eigenfunctions are 
broadened versions of the local solutions but their qualitative similarity is evident 
in Fig. H Such broadening of the eigenfunctions may be achieved from an anoma- 
lous, enhanced diffusion. For A in the range applicable to spiral galaxies, namely 
0.01 < A < 0.1, the eigenvalues F (growth rate) resulting from the non-local prob- 
lem are smaller, while the eigenfunctions are generally wider - see Fig. El As might 
be expected, the non-local effects become weaker as A decreases, i.e., as the disc 
becomes thinner. 




(37) 
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Figure 1: The first two eigenfunctions Qi and Q2 of the local and non-local dynamo 
equations (|T^ and (|TT|l respectively, for the case A = 0.01 and 7(r) of the form (jHZj) 
with 7m = 10, tq = 1 and R = 2. The local (non-local) eigenfunctions Q are shown 
with solid (dashed) curves and labelled with the superscript / (n). 
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-0.5 




—41 I \ I I I 

-7 -6 -5 -4 -3 -2 -1 

In A 



Figure 2: The dependencies of the leading eigensolutions on the dimensionless 
disc thickness A showing the variations on a log-log scale: (a) the width of the 
eigenfunctions Ar (i.e., the distance between positions where Q{r) decreases to e"^ 
of its maximum value) and (b) the decrement 7^ — T of the eigenvalues. The local 
and non-local solutions are distinguished by labels / and n, respectively. 
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o.z j-/iiiear suiuiiuiis lor iiie case > ^ 

When 7(r) is maximised off the disc axis, the eigenf unctions are locahsed in its 
vicinity on the length scale Ar oc A'^, while the corresponding eigenvalue decrement 
7m — r, which measures the growth rate F, scales as A^*. These two important 
qualitative measures of the non-local and local eigensolutions are distinguished by 
s = 1/3 and s = 1/2, respectively (see (jHUj) and Paper I). Though these scalings are 
derived for our present case > 0, they appear to hold also for the case = 
of Section 13.11 as Fig. |21 illustrates. Another important distinction of the non-local 
eigenf unctions is that they decay only algebraically at r oo (Soward, 1992a,b), 
in contrast to the exponential decay of the local eigenf unctions (e.g.. Moss et ai, 
1998). 

As an illustrative example we consider the local growth rate 



-fmr{2rm-r)/r^, r < 2rm , 
7(r) = { {3t 

7oo, r > 2rr, 



' m 



similar to (^Hj) with 7^ = — 27^/^^ giving e = (Ar^)-*^/^. The leading eigenfunction 
for the parameter values A = 0.01, 7^ = 1, 7oo = 0, = 0.25 and R = 4 is shown 
in Fig. El Other values of 700 are employed in Fig. IH below. 

We compare our eigenfunction with the asymptotic solution (jH^ valid close to 
on the length scale e = (Ar^)^/^, shown with dashed curve in Fig. El For 
e <^ |r — r^l <^ ^'m it determines the power law behaviour Q oc \r — rml'"^ (see (jSlj)) 
and this remains throughout the range < r < over which 7(r) has a quadratic 
profile. This actually links to a line dipole source, for which the poloidal magnetic 
field decays as |r — rm\~^- For r > 2rm, where 7(r) = 700 is a constant, we have 
instead (j2Sl), which takes the form 

XQ d'^G , , 

q{r) ^ rmTT^ — [r,rm) , (39) 

7m - 7oo drdrm 

where r^d'^G / drdrmir^rm) has the asymptotic representation (P7|) giving 

3Ar2 Q 

^ , " , , (40) 
^ 2(7m - 7oo)r4 V ^ 
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-5 -4 -3 -2 -1 1 2 3 4 5 

log2 (r/r^) 

Figure 3: The leading eigenfunctions Q{r) appropriate to 7(r) defined by with 
the parameter values A = 0.01, 7^ = 1, 7oo = 0, = 0.25 and i? = 4 for the 'local' 
p4|l (dash-dotted) and the non-local (fTT|) (solid) problems; remember that 7 = 700 
on 1 < log2 (r/r^) < 4. The local solution has an exponential tail for large \r — r^l, 
whereas the non-local solution exhibits an algebraic tail. The dashed curve shows 
the WKB asymptotic solution ()32|) with the scaling (jHOI). The dotted line shows the 
far field asymptote (fH^ or pUjl . valid for r ^ 2rm ■ 
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similar to ifzi'ji . iiiis is tne reason lor tne r power law eviaent m r ig. li^ior r > zr^ 
and corresponds to the algebraic decay of its quadrupole source. This asymptotics 
improves as joo decreases, and the exponent —4 was generally found in the numerical 
results independent of the value of 700, at least for 700 < 0. 

The main points to note are that the asymptotic solution generated by the 
Airy function approximates the numerical solution in the vicinity of r^- The nature 
of Fig. ini (namely, its relatively narrow relevant range in r) makes it impossible to 
easily identify the power law |r — r^l"^, given by (jSH), inside the source region 
< r < 2rm- Whether or not it is there, the power law outside the source 
region r > has a different explanation for its similar power law exponent of 
—4 given by (^Dj), specifically stemming from the quadrupole nature of the source. 
We also note that the horizontal and vertical components of magnetic field scale 
differently with r in the far field: B oc Q oc r~^, = —dA/dz oc r~^, and B^ = 
r~^d{rA)/dr oc dQ/dr oc r~^. Thus, the magnetic field becomes more horizontal 
with r at large distances from the dynamo active region. 

3.3 Steady states 

A simple heuristic form of nonlinearity in the radial dynamo equation (|T4|l . resulting 
from the standard from of a-quenching, a = ao (l + B"^ / Bq) ^, with ao the back- 
ground unquenched value, was derived by Poezd et al. (1993) who suggest that 7(r) 
is to be replaced by 

7(°)(r,t) = 7(r)[l-Q^(r,t)/g^(^)], (41) 

where Qoir) is essentially equal to Bq, but with the possible 2;-dependence of Bq 
averaged out. In this section, we compare the local and non-local steady states 
using this form with Qo{r) = 1 in both (|THl and (|TT|l . together with 7(r) of the form 
()38|1 where 700 was varied to explore properties of the solution at r > ri. 

Taking the eigenfunctions of Sect. 13.21 as initial conditions, we solved the Cauchy 
problem for both the local and non-local dynamo equations using the Runge-Kutta 
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Figure 4: Steady-state nonlinear non-local solutions obtained for 7(r) of the form 
with R = 4, = 0.25, A = 0.01, and = 0, -0.1, -0.3, -1.0 (with the 
value of 7oo shown near the corresponding curve). The algebraic tail remains similar 
to that in the linear solutions, Q ~ r~^, with the exponent independent of 700, 
although the magnitude of Q in the tail reduces as 700 decreases. 

method. This was performed also for the corresponding linear equation in order to 
verify the accuracy of the computations. 

Figure m shows steady states for the non-linear regime. The non-linear solution 
clearly has a broader distribution in r than the associated eigenfunction in Fig. Ob- 
The algebraic tail Q ~ characteristic of the non-local eigenmodes persists in 
the nonlinear solutions. Its functional form is independent of the local value of 7, 
as shown by our solutions with various values of 700. This is consistent with the 
fact that the algebraic tail results from the non-local magnetic connection with the 
dynamo- active region at smaller r. However, the magnitude of Q{r) at large r is 
affected by the value of 700. This effect can be estimated from with F = 0, 
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wnicn yieias q[r)(x y/7oo- J-M values oi q over tne appropriate range oi r m r ig. (j4j) 
agree with this to within 30% for = —0.3, and 20% for = —1. 

3.4 Propagating Magnetic Fronts 

If the local growth rate, 7o('"), is localized in a limited radial domain, growing 
magnetic field distribution of the local kinematic dynamo problem, (jlH), propagates 
at a finite velocity, V, into the dynamo-inactive region. This phenomenon was 
discussed by Moss et al. (1998, 2000), Petrov et al. (2001), Petrov (2002) and 
Fedotov et al. (2003). The eigenf unctions of ()14|1 decay exponentially at large r, 
Q{r) (X exp(rt — r/ro), where tq = (P/TY^"^ in dimensional variables. The position 
Tf where Q = const changes with time as rf = rgFt = (pvy^'^t, i.e., the magnetic 
front propagates at a constant speed V = {VPY^'^ (Moss et al, 1998). The front 
propagation is strongly affected by the asymptotic behaviour of Q{r) at r ^ 1: 
for the non-local solution, we have Q oc (r/ro)~^ exp (Ft), so that Q = const now 
yields rf = roexpFt/4, i.e., the magnetic front propagation is exponentially fast, 
with a speed increasing with time as oc expFt/4. This is a direct consequence 
of the non-local coupling of the inner regions, where the dynamo is strong, with 
the outer regions where the front propagates: since magnetic fields propagates at 
infinite speed through vacuum, magnetic field far from the dynamo region responds 
instantaneously to the exponential growth of magnetic field in the inner parts of the 
disc. 

These heuristic arguments are confirmed by numerical solutions illustrated in 
Fig. El where we show the position of the magnetic front defined here as a position 
r = Tf where Q{r{) = 10~^. The eigenfunction of the nonlocal problem was used as 
the initial condition, normalised such that its maximum was 10"^. The eigenfunction 
decreases with r relatively quickly for rm/2 < r < (within the dynamo region) 
and then slower at r > (where 7(r) = 700 < 0) — see FiglHl Correspondingly, the 
front first propagates slowly through the region where the solution rapidly decreases 
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1.5 




Figure 5: The position of a propagating magnetic front versus time for the hnear 
(sohd) and nonhnear (dashed) solutions of (fTT|) with 7(r) of the form (jHHj) with the 
same parameters as in Fig. and F = 0.852. The dotted line represents exp(Ft/4). 
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away rrom its maximum, out tne rront speeas up as it reacnes tne region or tne 
algebraic tail. Front propagation is inhibited by nonlinear effects as they suppress 
the exponential growth of the solution and then estabhsh a steady state. 

4 Conclusions and Discussion 

Arguably the most important property of the non-local solutions discussed above is 
the slow, algebraic decay of both kinematic and steady-state mean magnetic fields at 
large distances, along radius, from the field maximum. This feature can be important 
for the understanding of large-scale magnetic fields observed at large distances from 
the centres of spiral galaxies where both low gas density and reduced intensity of 
interstellar turbulence would preclude any strong magnetic fields. However, Han 
et al. (1998) have found indications, albeit not very strong, that the large-scale 
magnetic field in the disc of the Andromeda nebula (M31) at distances as large 
as 25kpc is almost as strong as at a radius of lOkpc. Beck (2003) argues that 
magnetic field in the galaxy NGC 6946 (one of the better studied nearby spiral 
galaxies) decreases with galactocentric radius unexpectedly slowly, so it can become 
dynamically important in the outer Galaxy. If approximated by an exponential 
profile, the energy density of the total (regular plus random) magnetic field has the 
radial length scale of perhaps as large as 8kpc. Meanwhile, the gas density decreases 
exponentially at a length scale of 2-3 kpc; if the turbulent velocity does not vary 
much with radius, this length scale applies to the turbulent energy density as well. 
Thus, the observations indicate that the large-scale magnetic field may decay with 
r slower than the turbulent energy density. 

The algebraic tail produced by the non-local magnetic efi'ects can be responsible 
for the slow decrease of the observed magnetic field with radius. It would be im- 
portant to test this idea with numerical solutions for a disc dynamo embedded in a 
poorly conducting halo. 

The assumption that the disc is surrounded by vacuum is essential for the nature 
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or tne long-aistance oenaviour or tne aynamo-generatea magnetic neia. tieai gaiactic 
discs are surrounded by turbulent halos, and the vacuum boundary conditions are 
justified if the time scale of the (turbulent) magnetic field diffusion above the disc 
is significantly shorter that through the disc, i.e., if 

Td /?h ' 

where /3h and /3d are the turbulent magnetic diffusivities in the halo and the disc, 
respectively, R is the radial length scale of the magnetic field in the disc, and L 
is the length of a magnetic line that leaves the disc through its surface and then 
returns to the disc with the radial distance R between its footpoints. As argued by 
Poezd et al. (1993), /9h//9d — 10-30, assuming that both the turbulent velocity and 
scale are in the halo are 3-5 times larger than in the disc. With L/R oi order unity, 
this yields Th/rd <^ 1; for L/R ~ 2 we obtain Th/ra ~ 0.1-0.5, a value arguably 
small enough to believe that the magnetic coupling through the halo is important. 

Another condition for the algebraic tail to be established is that the magnetic 
propagation time through the halo is short enough in comparison with the galactic 
lifetime, Th < 10^° yr. Using for L the radial length scale of a kinematic non-local 
magnetic mode in the disc, L ~ X''^l^h, i.e., assuming that the length of magnetic 
lines joining two positions in the disc through the halo docs not differ strongly 
from the radial separation of the positions), where h ~ 500 pc is the scale height of 
the galactic ionized disc and A = 0.05-0.1 is its aspect ratio, we obtain Th ~ (5- 
40) X 10^ yr, with /3h = (10-30)/3d and (3^ ~ 10^^ cm^ s~^ Considering the uncertainty 
of the halo parameters, this estimate seems to indicate that the algebraic tail can 
be estabhshed in spiral galaxies, especially in those with strong halo turbulence or 
small disc. 
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Appendix A. The evaluation of the eUiptic integral 
(1171) and its derivative 

To evaluate the elliptic integral ()17|1 . where we have introduced k via m = sin^ k, 
and used an iterative method described by Abramowitz and Stegun (1965): for the 
staring values 



ao = 1 



6o = cos K = Vl — m , 



(Al) 



we calculate 




'n+l — 




'n 



(A2) 
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to oDtam 



K(m) = £- . (A3) 



"oo 



The derivative dK/dm is calculated using a generalization of this algorithm sug- 
gested by W. Dobler (unpublished). We perturb the system above, m — >• m + 
5m, ttn ^ ttn — Kn5m, bn ^ bn — (n^m, whcrc 6m <^ 1. The starting values and 
Co are obtained from the Taylor expansions of (jAljl : 

= , Co 



and similarly from ()A2jl : 

1/ \ /• \ /• '^nCn ~l~ ^nb-n 

Perturbing ()A3|) with respect to m, we obtain 

dK vr 



*'oo 



dm 2a^ 

Typically, one or two iterations are sufficient to calculate K{m), dK/dm to the 
eighth decimal place. 



Appendix B. The regularization of the singular in- 
tegral kernel W{r,r'), (1151 ) 



Equation is solved numerically for < r < i?, with suitably large R. The inte- 
gral kernel of this equation, W{r, r') has a singularity along the line r = r'. Simple 
quadrature methods may display poor convergence if such singularities are ignored 
and do not apply at r = r'. These difficulties are resolved using the subtraction of 
the singularity (Press et al., 1993) by performing the following transformation: 

Ar/(r) 



[r-7(r)]g(r) 



W{r, r')Lr'{q{r'))dr' 



\r]{r) 



R 



W{r,r') Lriqir')) - Lriqir)) dr' + Lr{q{r))M{r 



(Bl) 
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wnere 



LM{r)) ^ 



d_ 

Or 



ld_ 
r dr 



M{r) 



R rR 

W{r^ r') dr' = r r' 
Jo 



Ji{kr)Ji{kr') dk 



dr'. 



(B2) 



As a result, the singularity in the integral operator at r = r' is regularized since 
Lr'{q{r')) — Lr{q{r)) = at r = r'. The additional term Xri{r)M{r)r~^Lr{q{r)) 
describes a modified radial diffusion with diffusivity Xri{r)M{r)r~^. 

For finite, non-vanishing values of r, M(r) can be calculated as follows. Making 
the substitution a = kr, we obtain 

rR 



M{r) 

and with x = r' /r we find 



rR/r 

M{r) = X 
Jo 



Ji{o')Ji{ar'/r) da 



dr' 



Ji(<7)Ji(crx) da 



dx = r'^M{r) 



where 



M(r) 



R/r 



W{l,x) dx 



is a universal function of r that need only be calculated once. For a mesh point r„, 
n < N, rN = R, we have 

"1 



M 



N 



W{l,x) dx . 







^ , rR/rn-i 

Mn-1 = Mn+ W{l,x)dx, 
jR/r„ 



n = N,...2 



(B3) 



where M„ = M(r„). Using an (open) quadrature method, we find M^- ~ 0.31409623, 
the remaining M„, n < N depend on the descretisation. From the recurrence relation 
()B3|1 . M„, and hence M„, may be precomputed by a single integration over the range 
1 to R/ri. 

We now take centered differences on the regularized equation ()B1|) . The point 
j = i in the corresponding sum is omitted since there [Lrj{q{rj)) — Lr^(g(rj))] = 0. 
This leads to the eigensystem FQ = AQ which we solve using the QR algorithm. 
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